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ABSTRACT 

Optical spectroscopy and photometry of SN 2006aj have been performed with the Subaru telescope 
at t > 200 days after GRB060218, the X-ray Flash with which it was associated. Strong nebular 
emission-lines with an expansion velocity of v ~ 7,300 km s -1 were detected. The peaked but 
relatively broad [OIJAA6300, 6363 suggests the existence of ~ 2M Q of materials in which ~ 1.3M Q 
is oxygen. The core might be produced by a mildly asymmetric explosion. The spectra are unique 
among SNe Ic in (1) the absence of [Call] AA7291, 7324 emission, and (2) a strong emission feature 
at ~ 7400A, which requires ~ 0.05M Q of newly-synthesized 58 Ni. Such a large amount of stable 
neutron-rich Ni strongly indicates the formation of a neutron star. The progenitor and the explosion 
energy are constrained to 18M < M ms < 22M Q and E ~ (1 - 3) x 10 51 erg, respectively. 

Subject headings: gamma-rays: bursts - supernovae: general - supernovae: individual (SN 2006aj) - 
nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

The connection between long soft Gamma-Ray Bursts 
(GRBs) and Type Ic Supernovae (SNe Ic) has been ob- 
servationally established (see Woosley & Bloom 2006 and 
references therein). It has been argued that well exam- 
ined GRB-associated SNe (GRB-SNe) are very energetic 
explosions, called hypernovae, resulting from a massive 
stellar death. The kinetic energy (Ek) is E 5 i = i?K/10 51 
erg s" 1 > 10, and the main-sequence mass of the pro- 
genitor is M ms ~ 4OM (e.g., SN 1998bw: Iwamoto et 
al. 1998). 

Recently, a new connection between the X-Ray Flash 
(XRF), which is a low energy analog of a GRB, and an 
SN Ic has been found. In association with XRF 060218 
(Campana et al. 2006), SN Ic 2006aj was discovered 
(Ferrero et al. 2006; Mirabal et al. 2006; Modjaz et al. 
2006; Pian et al. 2006; Soderberg et al. 2006; Sollerman 
et al. 2006). 

At early phase, SN 2006aj showed velocities interme- 
diate between GRB-SNe and canonical SNe Ic. The 
OIA7774 absorption was weak. The light curve evolved 
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more quickly than other GRB-SNe. Mazzali et al. (2006) 
concluded that SN 2006aj is less energetic (E51 ~ 2) 
than other GRB-SNe and ejected a smaller amount of 
56 Ni [M( 56 Ni) ~ 0.2M o ]. They suggested that the main- 
sequence mass of the progenitor is M ms ~ 20 — 25M Q , 
and hence that XRF-associated SN 2006aj was driven by 
a neutron star (NS) formation. 

In this Letter, we present observations of SN Ic 2006aj 
at t > 200 days with the Subaru telescope (hereafter 
t is the rest frame time since 2006 February 18). We 
show that the nebular phase data provide constraints on 
nucleosynthesis and the compact remnant. 

2. OBSERVATIONS AND DATA REDUCTION 

We performed spectroscopic and photometric observa- 
tions of SN 2006aj on 2006 September 17 (UT: t = 204 
days) and November 27 (273 days) with the 8.2 m Sub- 
aru telescope equipped with the Faint Object Camera 
and Spectrograph (FOCAS: Kashikawa et al. 2002). For 
spectroscopy, we used 0'.'8 width slit and the B300 grism, 
which gave a wavelength coverage of 4700-9000 A and a 
spectral resolution of ~ 10.7 A. The exposure times were 
2,800 s and 5,400 s, respectively. For flux calibration, 
GD71 and BD +28°4211 were observed. 

For photometry, we obtained B- (60 s), V- (60 s) and 
i?-band (60 s) images in September and B- (300 s), V- 
(180 s), and i?-band (180 s) images in November. We 
obtained images of standard stars around PG0220+132 
and SA98-634 for photometric calibrations. Our photom- 
etry is shown in Table 1 and the photometry-calibrated 
spectra in Figure 1. 

The spectra are the sum of the SN and host galaxy 
spectra (Fig. 1). We find broad emission lines with the 
corresponding expansion velocity v ~ 7, 300 km s" 1 at 
both epochs (Fig. 2: see also Foley et al. 2006). Espe- 
cially strong is [OIJAA6300, 6363. So the existence of a 
SN Ic in XRF 060218 is unambiguously confirmed. 

The continuum level is almost identical in both epochs. 
Since no strong emission from an SN is expected at ~ 
6100, 6800, and 8000A, we regard the flux level of the 
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continuum connecting these wavelengths to be the host 
contribution, expressed as —2.0 x 1CT 17 x (A/6000A) + 
4.1 x 10~ 17 erg s _1 cm~ 2 A -1 (Fig. 1). The derived host 
magnitude is fainter than that estimated by Sollerman et 
al. (2006) by about 0.3 mag in V- and i?-bands, while the 
color V — R is consistent. The difference probably stems 
from the systematic errors. The host is blue, as expected 
from its faintness and low metallicity (Sollerman et al. 
2006). 

3. ANALYSES OF SPECTRUM FEATURES 

Figure 2 shows that the nebular spectra of SN 2006aj 
have the following two unique features compared with 
other SNe Ic. (1) [CaII]AA7291, 7324 is absent. (2) A 
strong emission feature is seen at ~ 7, 400A. In this Let- 
ter, we adopt ^ = 35.84 and E(B — V) = 0.13 (Pian 
et al. 2006). Our discussion is based on the September 
spectrum. 

3.1. Total Mass at v < 7, 300 km s" 1 

The bolometric luminosity (Lboi) is expressed as (e.g., 
Maeda et al. 2003a), 

iboi = (D 1 + 0.035)L 7 , (1) 

at a nebular epoch. Here the 7-ray energy input rate 
from M( 56 Ni) = 0.2M© is L 7 = 4.4 x 10 41 erg s" 1 at t ~ 
204 days. Z? 7 (<~ 3r 7 /4 for a homogeneous nebula, which 
is a good approximation at low velocities) is the fraction 
of the 7-ray energy deposited within the SN nebula, and 
0.035 accounts for the in situ positron deposition rate. 
The optical depth to the 7-rays is 

Ti = 0.071 ( 7j300kms -i) (^04d^i) ' 

(2) 

where M7300 denotes the total mass ejected with v < 
7,300 km s" 1 . Substituting L hol ~ 6.2 x 10 40 erg s" 1 
(taking into account the reddening and a typical NIR 
contribution ~ 20%; Tomita ct al. 2006) into equations 
(1) and (2), we obtain M 730 o ~ 2M . 

3.2. Absence of Ca Lines and Ejecta Structure 

No detection of [Call] and Call IR triplet provides the 
following constraints. In a core-collapse SN explosion, 
more than 4 x 1O~ 3 M of Ca is synthesized if -E51 > 1 
(e.g., Nakamura et al. 2001). In an SN Ic, the emissivity 
of [CaII]AA7291, 7324 is about three orders of magnitudes 
larger than [OI]AA6300, 6363. If Ca co-existed with O, 
Ca + would be the predominant ion (e.g., Fransson & 
Chevalier 1989). Thus, if Ca (> 4 x 10~ 3 M Q ) and O 
were microscopically mixed, then Mo would have been 
much larger than 20M Q to satisfy the observed upper 
limit for the [Call] (Fig. 2). This is inconsistent with 
M 730 o ~ 2M (§3.1). 

The separation of Ca and O occurs naturally in or- 
dinary supernova models. This is also seen in Spitzer 
observations of Cas A (Ennis et al. 2006). Nucleosynthe- 
sis calculations show that the explosion produces chem- 
ically distinct regions: Fe-rich materials (where 56 Ni is 
the predominant product), Si-rich materials (where most 
of Ca is produced), and O-rich materials (mostly O and 
some unburned C). These regions may well be mixed 



macroscopically in velocity space (e.g., by jets: Maeda & 
Nomoto 2003b; Maeda 2006b), but not microscopically. 
Hereafter, M Fe _ ri ch, M S i_ r i c h, and M - r ich, respectively, 
denote the masses of the Fe-rich, Si-rich, and O-rich ma- 
terials ejected with v < 7, 300 km s _1 . 

In this "microscopically separated" model, the sum of 
the [Call] and Call luminosities (Lcaii) is related to the 
mass of the Si-rich region rather than the mass of Ca 
itself, i.e., 

iCall ~ ^Asi-rich-L^asi/call , (3) 

where the fraction of the absorbed 7-ray energy go- 
ing into the Ca lines would be /call ~ 0.5, with the 
remaining fraction mostly going into [Sil] NIR lines. 
Asi = Msi-rich/M7300 is the mass fraction of the Si-rich 
region. A factor, as;, accounts for the fact that the 7- 
rays are scattered more easily within the Fe-rich region 
than in the other regions, because the 7-rays begin their 
flight in the Fe-rich region. At late epochs, the mean free 
path of the 7-rays is larger than the size of the nebula, 
thus asi is an order of unity. Taking the upper limit 
Lean < 7 x 10 39 erg s _1 , we obtain 

M Si _ rich <O.6M (l/asi)(O.5// Ca ii) . (4) 

The upper limit is even tighter according to the spectrum 
synthesis (§3.5), i.e., Msi- r i c h < O.15M , because other 
lines (e.g., [Fell]) should also have large contribution. 

3.3. Masses of Oxygen and Fe-peak Elements 

From the masses we have obtained as M7300 ~ 2M , 
^Si-rich < O.15M , and M( 56 Ni) - 0.2M Q , we derive 
Mo-rich ~ 1.6M© and M Fc _ ri ch ~ 0.3M Q . We here 
examine if Mo -rich and Mpc-rich are consistent with the 
reddening-corrected luminosity of [OI]AA6300,6363 (Lqi) 
and that of [Fell] and [Nill] (L Fo -pcak: integrated in 4700 
- 5500A and 7000 - 7600A). 

Following the argument similar to the one used in de- 
riving equation (4), we obtain 

Mo- r ich ~ 1.6M (O.7/a o )(O.7//oi) , (5) 

with L i ~ 1-7 x 10 40 erg s _1 . Hereafter, a; and A; 
(where i = Fe, Si and O) are used for the Fe-, Si-, and 
O-rich regions in the same way as asi and Asi- r ; c h- 
/oi is the fraction of the absorbed 7-ray energy go- 
ing into [OI]AA6300,6363 (with remaining fraction going 
into [CI]A8727 and OIA7774). Equation (5) shows that 
ao ~ 0.7 is necessary to reproduce Loi- To meet the 
trivial relations A;Q!i = 1 and A; = 1, a-p c ~ 2.6 is 
required. With this value, the spectrum synthesis (§3.5) 
well reproduces LFc-pcak ~ 2 x 10 40 ergs s _1 . 

These values for a; are physically reasonable. For an 
SN nebula having a constant density with MF C _ r ich = 
O.3M , M Si _ ri ch = 0.1M Q , and M _ r ich = 1.6M , we 
derive a Fc -rich ~ 2.9, a si _ ri ch — 1.8, and a -rich ~ 0.6 
by our ID 7-ray transport calculation. These are ob- 
tained for a concentric distribution (with the Fe region 
in the center), but it gives an approximately correct es- 
timate even for a jet/blob-likc distribution of the Fe-rich 
materials. 

3.4. [NUIJX7380 and stable 58 Ni 

The strong emission at ~ 7400A is not common for 
SNe Ib/c, although in other SNe Ib/c the feature may be 
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hidden by the [Call]. The contribution from [FeII]A7388 
and [FeII]A7452 (a 4 F-a 2 G) should be small. The radia- 
tive transition probabilities of the [Fell] are smaller than 
those of [FeII]A7172 and [FeII]A7155 that arise from the 
same transitions, and the luminosity at 7100 — 7200A is 
only one-third of the luminosity of the 7400A feature. 

We suggest that the 7400A feature is [NiII]A7380 (Fig. 
2), as also seen in the Crab Nebula (MacAlpine et al. 
2007). From the spectrum synthesis (§3.5), we find that 
this feature is reproduced by [NiII]A7380 if the mass of 
Ni is ~ 0.05M©. Since virtually all 56 Ni already decayed 
at t > 200 days, this must be stable 58 Ni. 

3.5. Spectrum Synthesis 

To confirm the above results, we perform nebular 
spectrum synthesis calculations. We set the masses at 
v < 7,300 km s _1 as follows: Mo- r ich = 1.6M (where 
M = 1.3M© and M c = O.3M ), M si _ r ich = O.1M 
(where M si = 0.08M©, M s = 0.01M©, M Ca = 4 x 
1O~ 3 M , and M Fo = 4 x 1O~ 3 M ), and M Fc _ ric h = 
O.3M [where M( 56 Ni) = O.2M and M( 58 Ni) = 
O.O6M ]. We set a Fc = 2.6, a Ca = 1-0, and a = 0.7 for 
the 7-ray deposition (§3.3). 

A one-zone spectrum synthesis is applied to the O- 
rich, Si-rich, and Fe-rich materials. A Monte-Carlo 7-ray 
transport calculation is performed for a uniform nebula 
with Af 730 o = 2M , M( 56 Ni) = O.2M . The deposited 
luminosity is apportioned to each region according to 
X\a\ (equation 3). Positrons deposit their energy only 
to the Fe-rich materials. 

With the deposition luminosity, non-LTE rate equa- 
tions coupled with the ionization-recombination equilib- 
rium are solved for each region (Mazzali et al. 2001). The 
model spectrum fits the observed one fairly well (Fig. 3). 
The feature at ~ 7, 400 A is well explained by [NiII]A7380 
(§3.3) from - O.O5M of stable Ni. 

4. PROPERTIES OF THE PROGENITOR AND ITS 
EXPLOSION 

We have derived the masses as follows (§3): M7300 ~ 
2M , Mo-rich ~ 1.6M (i.e., M ~ 1.3M if 80% is 
oxygen), M S i_ r ich < O.15M , and M Fc _ ric h ~ O.3M 
(where M( 56 Ni) - O.2M and M( 58 Ni) - O.O5M ). 

Progenitor Mass: The masses of Mo ~ 1.3M and 
-W7300 ~ 2M provide strong constraints on M ms as fol- 
lows. Theoretical models show Mo < 1.3M in the 
whole ejecta if M ms < 18M (Nomoto & Hashimoto 
1988). This leads to the lower limit of M ms > 18M . 

For given £51, larger M ms results in larger M^oo- Our 
spherical calculations show M7300 > 2M if M ms > 
22M and -E51 = 2. Since explosion asymmetry only in- 
creases the mass at the low velocity (Maeda et al. 2002; 
Maeda & Nomoto 2003b) , the result of the spherical cal- 
culations can be used to place the upper limit to M ms . 
We thus obtain M ms < 22M . Therefore, we narrow 
down the progenitor mass range as 18M < M ms < 
22M . 

The total oxygen mass is ~ 2.3M in the ejecta, since 
the oxygen mass at v > 7, 300km s _1 is ~ 1M (Mazzali 
et al. 2006). This is consistent with M ms ~ 22M . 

Explosion Energy: M Si _ rich < O.15M (§3.2) gives 
information on E51. The mass of the Si-rich region is 



larger for larger M ms (because of the higher progeni- 
tor density) and for larger E$i (because of the more 
extended burning). According to our explosion models 
(e.g., Nakamura et al. 2001), we obtain the constraints 
of M ms < 25M and E 51 < 3. 

NS Formation: M( 58 Ni) ~ O.O5M provides a 
strong indication of the NS formation. The large amount 
of 58 Ni indicates that neutron-rich materials are ejected 
from the vicinity of the central remnant. Our models 
for M ms = 2OM and £51 - 2.5 shows that M( 56 Ni) 
~ O.2M and M( 58 Ni) ~ O.O4M if the remnant's mass 
is M cut ~ 1.4M . If M cut ~ 1.5M , then M( 56 Ni) ~ 
O.16M andM( 58 Ni) - 9xlO~ 4 M . For M cut ~ 1.6M , 
M( 56 Ni) - O.O7M and M( 58 Ni) - 3 x 1O- 4 M . The 
NS formation is also consistent with M( 56 Ni) ~ O.2M . 

5. DISCUSSION 
5.1. Signatures of an Asymmetric Explosion 

In GRB-SN 1998bw, the observed mass at low veloc- 
ities exceeds that predicted by the spherical explosion 
model, which suggests that the explosion is aspherical 
(Mazzali et al. 2001; Maeda et al. 2003a, 2006ac). On 
the other hand, the existence of the asymmetry is not 
clear in SN 2006aj. M 730 o ~ 2M (§3) is consistent 
with the spherical explosion model for a CO star with 
M ms ~ 22M and E 51 ~ 2. 

We note, however, that the spherical model has a low 
density hole at v < 2, 000 - 2, 500 km s _1 . The mildly- 
peaked [OI]AA6300,6363 profile shows an enhancement of 
the material density at v < 3, 000 km s . This indicates 
an asymmetric explosion, as also suggested by polariza- 
tions (Gorosabel, et al. 2006). The asymmetry might be 
produced either by a jet-like explosion or Raylcigh- Taylor 
instability at the Si/CO interface (Kifonidis et al. 2003). 
If this is a signature of a jet-like aspherical SN explosion, 
the jet would be wider than in SN 1998bw (intrinsically 
or due to stronger lateral expansion; Maeda & Nomoto 
2003b), since the signature is seen only in the innermost 
part (see also Mazzali et al. 2007). 

5.2. Relations to other SNe Ib/c and SN 1998bw 

The absence of [Call] AA7291, 7324 indicates that the 
chemically distinct regions are not microscopically mixed 
to one another (§3.2). This is probably the case in SNe 
Ib/c in general: otherwise the [Call] would be much 
stronger than [OI]AA6300, 6363 in most SNe Ib/c, as is 
contrary to observations (e.g., Matheson et al. 2001). 

Why the [Call] is absent only in SN 2006aj may be 
answered from a tendency of E^i as a function of M ms . 
For progenitors with M ms < 2OM , the [Call] can be 
strong relative to the [OI], because Mo-rich is smaller 
than in SN 2006aj despite similar Msi_ r ich (since E^i ~ 
1-2). For hypernovae with M ms ~ 4OM (e.g., SN 
1998bw), the [Call] can be strong since the large £51 > 
10 yields large Msi- r ich- Thus, presence of the [Call] 
in the majority of normal SNe Ib/c suggests that such 
normal SNe Ib/c may have originated in relatively low- 
mass stars with M ms < 2OM . This implies that they 
evolve through binary paths, since single stars with this 
mass range are expected to end as SNe II (Hirschi et al. 
2005). 
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TABLE 1 
Photometry 



Date/ Component 


B 


V 


R 


September 17 a 


20.64 ±0.1 


20.23 ±0.1 


19.82 ±0.1 


November 27 a 


20.71 ±0.1 


20.34 ±0.1 


19.86 ±0.1 


Host (Subaru) b 




20.53 


20.16 


SN on September 




22.45 


21.74 


SN on November 




22.94 


22.21 


Host (Sollerman) d 


20.46 ± 0.07 


20.19 ± 0.04 


19.86 ±0.03 



a Using the FOCAS filter system, which is similar to the Johnson-Cousins system (Kashikawa ct al. 2002). 
^See §2 and Figure 1 for the assumed host spectrum. The photometry is given for the Johnson-Cousins system. 

c Thc intrinsic SN magnitude obtained after subtracting the host contribution and removing narrow emission lines. In the Johnson-Cousins system. 
^The host magnitude estimated by Sollcrman et al. (2006) 
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Fig. 1.— Reduced spectra of SN 2006aj on 2006 September 17 (red), and on 2006 November 27 (blue). The flux is calibrated with the 
photometry. The redshift of the host (z = 0.0335) is corrected for. The assumed host galaxy contribution (§2) is shown by the black dotted 
line. The intrinsic SN spectra, after subtracting the host contribution and removing major narrow emission lines, are also shown (bottom). 
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Fig. 2.— Comparison of emission profiles, (a) Comparison of [OI]AA6300, 6363 (black) with [Call] AA7291, 7324 (red). The center of the 
wavelength is assumed to be 6300 and 7308A for the [OI] and [Call], respectively, (b) Comparison of the [OI] (black) with the emission 
feature centered at 7380A(blue). We suggest that the feature is [NiII]A7380. Note that an atmospheric absorption at ~ 7, 400A is not 
removed. 
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Fig. 3. — Synthetic spectrum as compared with the September spectrum. The model spectrum, reddened with E(B — V) = 0.13, is 
compared with the host-subtracted SN spectrum at the SN rest frame (blue). For presentation, the fit to the original spectrum is also 
shown (red), in which the host spectrum (Fig. 1) is added to the model spectrum. For the model spectrum, the individual contributions 
from the O-rich (black solid), Si-rich (dashed), and Fe-rich (dotted) regions are shown. 



